Relativistic treatment of spin-currents and spin-transfer torque 
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It is shown that a useful relativistic generalization of the conventional spin density s(r, t) for the 
case of moving electrons is the expectation value {T{f, t), T4,(f, t)) of the four-component Bargmann- 
Wigner polarization operator r,i = (T, T4,) with respect to the four components of the wave function. 
An exact equation of motion for this quantity is derived, using the one-particle Dirac equation, and 
the relativistic analogues of the non-relativistic concepts of spin-currents and spin-transfer torques 
are identified. In the classical limit the time evolution of (r, T4,), the integral of (T(r, t), 74(r, t)) over 
the volume of a wave packet, is governed by the equation of motion first proposed by Bargmann, 
Michel and Telegdi generalized to the case of inhomogeneous systems. In the non-relativistic limit 
it is found that the spin-current has an intrinsic Hall contribution and to order a spin-orbit 

coupling related torque appears in the equation of motion for s{f,t). The relevance of these results 
to the theory of the intrinsic spin Hall effect and current-induced switching are briefly discussed. 

PACS numbers: 71.70.Ej, 72.25.-b, 75.60.Jk, 85.75.-d 
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Central to the emerging technology of spin-based elec- 
tronics, often referred to as spintronics, is the observation 
that electron transport can be influenced not only by cou- 
pling to the charge but also to the spin of the current car- 
rying electrons, jlj A striking example of charge-current 
being effected by the magnetic state of the conductor is 
the giant magneto-resistance (GMR) phenomenon. Evi- 
dently, the complementary effects of charge-currents in- 
ducing changes in the magnetization of the conductors 
are also of interest. An example of these is the current- 
induced switching first predicted by J. C. Slonczewski 
and independently by L. Berger. [3, As is now well 
established, it is due to the spin-transfer torque that a 
spin-polarized current can exert on the magnetization of 
the structure through which it flows. P| However, the de- 
tails of how such torques come about have not been, as 
yet, fully explored. In particular, all discussions of the 
problem are currently based on non-relativistic quantum 
mechanics and hence neglect the spin-orbit coupling. The 
purpose of this letter is to present a fully relativistic the- 
ory of spin-currents and the above spin-transfer torque 
in order to provide a conceptual framework in which the 
spin and orbital degrees of freedoms can be treated on 
equal footing. 

Another topic, to which such developments are rele- 
vant, is the spin Hall effect intensively studied in semi- 
conductor spintronics. 0] It involves a spin-current 
flowing perpendicularly to a charge-current in a sample 
of finite width. Interestingly, it implies spin accumula- 
tion at the edges and the possibility of spin-injection into 
an adjacent sample without the presence of magnetic or 
exchange fields. Here, spin-orbit coupling is the central 
issue and a source of difficulty is the lack of a well-defined 
spin-current in a spin-orbit coupled system. 6] It is hoped 
that the polarization-current introduced in this letter will 
clarify this matter considerably. 



A moving electron carries with it a spin and this mov- 
ing spin amounts to a spin-current. Classically, it is given 
by the tensor product, J ci = Sd v, of the velocity vec- 
tor V and a classical spin vector Sd- Quantum mechani- 
cally, for an electron described by a two-component wave 



function (h = 



with spin components 0| and 



it is given by the tensor density 
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where J = ih{V — V) / (2mc), the wave functions and 
therefore also the spin-current J s are evaluated at the 
space time point (r, t). The physical significance of J s 
becomes apparent if we study the time evolution of the 
spin density defined by s = (jf^ S (f>. From the time- 
dependent Schrodinger equation which includes a Zee- 
man term of the form —j-iB^ ■ B we find 
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Clearly, V • J s may be regarded as a torque which, in ad- 
dition to the more familiar microscopic Landau-Lifshitz 
torque sxB, causes the spin density sat the point to 
evolve in time. As explained, at length, in the insightful 
review of Stiles and Miltat this spin-transfer torque de- 
pends linearly on the charge-current and plays a central 
role in current-induced switching. 

Note that for i? = Eq. ||2J) is a continuity equation for 
the spin density s and as such it follows via the Noether 
theorem from the fact that s is a conserved quantity. The 
difficulty of generalizing the continuity equation to spin- 
orbit coupled systems, such as described by the Dirac 
equation or by various model Hamiltonians, e.g. , like 
those used in semiconductor spintronics, 0j0] arises from 
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the circumstance that in these cases the spin operator no 
longer commutes with the Hamihonian and hence the 
conventional spin density is not conserved as the time 
evolves. In what follows this dilemma is resolved by the 
choice of a convenient, covariant description of the spin- 
polarization, as an alternative to that afforded by the 
usual spin operators in non-relativistic quantum mechan- 
ics. In the following, for the sake of simplicity only the 
case of non-interacting positive energy Dirac electrons 
will be considered. 

For electrons described by the Dirac equation in the 
standard representation it is common to refer to 
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as the 4x4 Pauli spin operator. However, it corre- 
sponds to the spin of an electron only in its rest frame 
and hence its use is not convenient in the case of many 
moving electrons. Moreover, as mentioned above, it does 
not commute even with the field-free Dirac Hamilto- 
nian — ca. ■ p + (3mcC^ and hence the correspond- 
ing density is not that of a conserved quantity. A more 
suitable approach for describing the spin-polarization of 
moving electrons is to use the four-component polariza- 
tion operator = (T^T^) introduced by Bargmann and 
Wigner. The salient features of this approach and its 
relations to other, alternatives, are fully discussed in a 
comprehensive review article by Fradkin and Good. 
Here they will be merely referred to as the need arises. 

For the case of one electron, in the presence of an 
electromagnetic field described by the vector potential 
A = A{r,t) and a scalar potential V — V{r,t), the four- 
component polarization operator is defined by 



f = /3S-iS4 — 
rricC 



(4) 



T4 = i S • 
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where the canonical momentum operator takes its usual 
form: II — {p — eA)l4, with 1 4 being the 4x4 unit ma- 
trix, S is the spin operator defined in Eq. For future 
reference note that S is part of the four-component op- 
erator = (S, E4) whose 4th component is defined as 
E4 = — i 75, see e.g. Ref. Iioi It is also of interest to note 
that both and are covariant axial four-vectors. 

From the point of view of our present concern the 
most important property of is that commutes with 
the field-free Dirac Hamiltonian. Thus, as will be shown 
below, the corresponding vector density satisfies a conti- 
nuity equation. To see the connection between the non- 
relativistic spin operator a and T it is useful to note that 
the latter is related to the magnetization of an electron 
in its rest frame by a Lorentz boost. 



To derive a relativistic analogue of Eq. ^ one has to 
calculate the first derivative with respect to the time of 
the polarization densities T — T{f,t) and T4 — T4{f,t) 
defined by 

f = tp+ftl; and T4 iIj+T4^P , 

where -0+ — ip^{f,t) is the adjoint (conjugate transpose) 
of the four-component solution -ijj = ipif^t) of the time- 
dependent Dirac equation corresponding to the Hamil- 
tonian Hu = Wd(^, — ca ■ U + (3mcC^ + eVl4. By 
using the chain rule for all four components fi = 1, . . . , 4 
in Eq. 

dt dt '^'^^'^ dt ^ ^ ^ dt ' 



and the relations 
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after some lengthy but straightforward algebra one ar- 
rives at 



V- J' + VJ" 



— Sx B 



■ES4 



and 



dT4 



rricC 



■S-E , 



(5) 
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where B = B{r,t) is the magnetic induction vector and 
E = E{r,t) the electric field intensity. In Eqs. lO and 
©, the four-component density = (5, S4) is given by 

5 = 54 = V'^E4V' 

and the polarization-current density tensors are defined 
as 

y' = cV'+(r®a)v J4' = c0+(T4a)V' 
J" = ^+(22/375)^- J a" 



0^+ 2 — X cE -0 



(7) 



such that 



Although the above relations bear some resemblance to 
Eq. (O the problem of generalizing it to relativistic quan- 
tum mechanics is not yet completed because Eqs. (jSJ and 
(0 ) are not a closed set of equations for the polarization 
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density in terms of corresponding currents. To proceed 
further one must derive the equation of motion for the 
four-component auxihary density iS^ . Following the same 
route as in the case of 7^ = (T, T4) one straightforwardly 
finds that 



find 



dS 



Tfi C — ' —> 

-icVSi = ^J^" + iV X J (8) 



and 



I — cV • 6 = « 

at h 



(9) 



with J — cip'^aip being the relativistic probability cur- 
rent density. Remarkably, some of the same currents 
which appear in the equations for 7^ determine 5^ and 
hence for a given set of currents in Eq. 0, Eqs. Q, 
© and lO can be solved for 7^ and iS^. 

The relations in Eqs. © - ® are the central result 
of this letter. Namely, the comparison of Eqs. l(2Jl and 
(O, without electromagnetic fields, uniquely identifies 
the polarization-current density a,s J ' + J"l3 and its 
divergence as the relativistic generalization of the con- 
ventional spin-current density and spin-transfer torque, 
respectively. Indeed, in the case of a vanishing electro- 
magnetic field Eq. ((S} reduces to a continuity equation 
for the polarization density in the same manner as Eq. ^ 
is a continuity equation for the magnetization density. 

To shed light on the physical content of these results 
they will now be examined in two separate limits. First 
the classical, ?i — > 0, then the non-relativistic, 
limit will be studied. 

The aim of the classical limit is to find a dynamical 
description of the polarization of an electron whose or- 
bital motion is classically given by the position vector 
fci{t) as prescribed by a relativistic classical mechani- 
cal equation of motion. This is the case of interest in 
Ref. 1^ and phenomenologically is treated in Ref. In 
short one assumes that ^/'(r, t) describes a wave packet 
centered at the position vector relit) and moving with a 
velocity Vd- The size of the wave packet is to be taken to 
be large compared with the Compton wave length h / mc 
but very much smaller then the scale on which the ex- 
ternal electromagnetic field vary. Moreover, it contains 
momentum components only near rricVci- In the lights 
of these assumptions it is natural to define a four-vector 
= (f, T4) as the integral of the density 7^ over the 
variable r within a volume fl as 



= f{t)= / d\f{r,t) , 



T4 



and derive an equation of motion for Tp(t) from those for 
7^ as given in Eqs. © and 10. Evidently, the property 
T^{t) is to be associated with a classical particle described 
by rc\(t) and v(t). Following the procedure of Ref. we 



j-^ f X B{f,i,t) (10) 



- i 7 ^ ^4-5(^:1, t) , 

TOcC 

where 7 is the Lorentz factor (l — w^/c^) and , as 
indicated, the external fields are evaluated at the cur- 
rent position rci(i) of the particle. Noting the result 
{Sfi) ~ 7~^T-^ of Ref. where (5^) denotes the clas- 
sical limit of the density 5^ integrated over fi, the form 
of the right-hand side of the above relation readily follows 
from Eq. jSJ. The term involving J^" was also shown to 
be zero by following the procedure of Ref. 9.. 

Before commenting on its most interesting feature, 
namely the torque on the left-hand side, it is reassur- 
ing to note that without it Eq. H10|) is exactly that of 
Bargmann, Michel and Telegdi (BMT) [l^ discussed at 
length by Landau. The extra factor of 7 is due to 
the fact that the derivative on the left-hand side is with 
respect to the global time t and not the pro per time in 
the rest frame of the electron as in Ref. [llj. To clarify 
the physical content of the BMT equation. Landau in- 
troduces a classical spin polarization vector in the rest 
frame, here denoted by s(t), and demands that when it 
is Lorentz boosted into the global frame, s{t) would be 
equal to f{t) in this letter. This allows him to convert the 
equation for f{t) to one describing s(t). In the interest 
of brevity his result is recalled here only to order c^^: 



ds 
'dt 
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Clearly, the second term on the right-hand side is a 

torque - due to an effective magnetic induction vector 

E X V I c - & moving electron experiences in an electric 

field E. In other words it is the simplest manifestation 

of spin-orbit coupling. For example, a Rashba Hamilto- 

nian analogue of this term is the main mechanism of the 

intrinsic spin Hall effect in the work of Sinova et al. 

The novel and the most interesting feature of Eq. Hl()|l 

is the torque on the left-hand side. For a uniform J ', 

< — > 

e.g. , in a bulk material, V • J7 ' is zero. But when the ef- 
fective classical particle crosses an interface or a domain 
wall with spin-dependent properties, the volume integral 
can be converted into an integral over a closed surface. 
If this surface includes the interface, whose differential 
oriented surface is dA, then the torque is ( — J '- 
) ■ dA, where ± refers to the opposite sites of the in- 
terface. Thus, just as in the non-relativistic case, see 
Ref. 5, this torque is due to an excess in the spin-current 
flowing in and out of the surface region causing a polar- 
ization accumulation. This then supports in very explicit 
terms, the identification of J ' + J " as the fully rel- 
ativistic generalization of the non-relativistic definition 
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of spin-current. Notably, in the classical limit one might 
approximate J ' ^ t ® u, where u is the relativistic ve- 
locity vector 7w and r is described by the BMT equation 
pO(l . This is a rather satisfactory result given the form of 
the non- relativistic classical formula J' ' s ® v quoted 
at the beginning of this letter. 

Finally, in what follows we examine the lowest order 
corrections to the non-relativistic theory. Working to the 
order 1 / c, after considerable algebra we find that the 
equations for T and S satisfy one and the same equation 
of motion, and the polarization-current is given by 



eA 
mo 



t/ (1) =0+ {S(E,J^ - 

^ + X I . (11) 



2mo 

Evidently, the first two terms in Eq. (|ll|l are the gener- 
alization of the conventional J s in Eq- O to the case of 
a non-vanishing vector potential. The third term 5Sf 
is a consequence of the internal contribution to the prob- 
ability current density due to the moving dipole moment 
d/int = hV X (0+0^0) / (2mo), 13 and its form readily 

follows from the anzatz Sj'^^^ = 4>^{B ® 6Jint)4>- [f^ 
Moreover, just as SJi^t does not contribute to the di- 
vergence in the continuity equation for the probability 
density, V • 6 J is identically zero and gives rise to no 
torque. 

To the order of 1 / there are many more contribu- 
tions. These will be discussed in a separate publication. 

^ — ^ (2) 

Here only one term 5 j7goc, which is clearly to be asso- 
ciated with the spin-orbit coupling, is highlighted: 



" soc 



2m'ic 



■E ■ 



(12) 















where (j) \s 4> renormalized as in Ref.[lll 

Remarkably, the off-diagonal terms have the form re- 
quired by the spin Hall effect, "i*] It means that for ex- 
ample an electric field, and presumably a charge-current, 
along the z axis, a spin polarization along the x axis 
implies a polarization-current in the y direction. In- 
terestingly, this term is the only contribution obtained, 
if one uses, in a simple minded derivation, the anoma- 
lous velocity, [l^ Wa — -eh{B x E) / {Amlc^), and 
the non-relativistic definition of the spin-current density 
® V + [v ® H Q Thus the first term in 

Eq. 11211 is a nontrivial consequence of our more general, 
fully relativistic treatment of the polarization. Clearly, it 
implies that for a charge current along the electric field 
there will be a helicity dependent contribution to the 
spin-current. 



Whilst the above reference to the spin Hall effect can 
not be taken as an explanation for the observed spin Hall 
currents due to the smallness of the vacuum coupling con- 
stant Asoc = eS/ (2mpC^), see Ref. 0, the presence of a 
relevant term in Eq. (|12|l suggests that an intrinsic spin 
Hall effect is a generic feature of spin-orbit coupled sys- 
tems and therefore of the relativistic quantum mechanics. 

In conclusion it should be stressed that the above dis- 
cussion was confined to a one-electron theory based on a 
one-electron Dirac equation. Nevertheless, the results es- 
tablish the line of reasoning a relativistic generalization 
of the corresponding many-particle theory has to take. 
In particular it will lead to a relativistic version of the 
semi-classical transport theory for the current-induced 
switching dynamics, 0] or for that of the spin Hall ef- 
fect. ^3 It will also facilitate the corresponding general- 
ization of the time-dependent density functional theory 
of Capelle et al., ^3 and will readily provide a frame- 
work for first-principles calculation usingftiUy relativistic 
methods such as the screened RKKR. (l^ Interestingly, 
it will also enter the relativistic Fermi liquid theory of 
Baym and Chin, [l^ 
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